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Abstract. We investigate the properties of the Gibbs states and thermodynamic 
observables of the spherical model in a random field. We show that on the low- 
temperature critical line the magnetization of the model is not a self-averaging 
observable, but it self-averages conditionally. We also show that an arbitrarily 
weak homogeneous boundary field dominates over fluctuations of the random field 
once the model transits into a ferromagnetic phase. As a result, a homogeneous 
boundary field restores the conventional self-averaging of thermodynamic observ- 
ables, like the magnetization and the susceptibility. We also investigate the ef- 
fective field created at the sites of the lattice by the random field, and show that 
at the critical temperature of the spherical model the effective field undergoes a 
transition into a phase with long-range correlations ^ r^"''. 



KEY WORDS: Critical fluctuations; disordered spin systems; Gibbs states; self- 
averaging. 



1 Introduction. 

The spherical model |5| is a lattice model where a (thermodynamic) random variable 
Xj is attached to every site j of a subset Vn of a (i-dimensional square lattice Z'^. 
This model is one of a handful of models where exact results can be obtained in 
the presence of a random field {hj,j G Z'^}. Thermodynamic properties of such a 
disordered spherical model outside the low-temperature critical line were studied by 
Pastur in the paper [13]. The magnetization on the critical line was also derived 
there in the limits — > ±0, where ho is the expected value of the random field. 

Some thermodynamic characteristics have discontinuities on the critical line, 
and, depending on the boundary conditions and the exact details of passing to the 
thermodynamic limit {Vn T ^'^); those characteristics can have different limiting 
values. Their values in the limits ho —>■ ±0 are, in some sense, extreme points of the 
sets of all possible limiting values. For some models those sets contain simply all 
linear combinations of the extreme values. For disordered models, like the spherical 
model in a random field, that is not necessarily the case. The aim of this paper 
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is to study thermodynamic properties of the spherical model directly on the low- 
temperature critical line. 

Many models in statistical mechanics are complicated enough to force us to 
restrict the investigation to finding only certain thermodynamic averages. For in- 
stance, sometimes investigation of magnetization is reduced to calculation of the 
averages 

jev 

where (■) denotes the average over the Gibbs distribution. However, as a rule, for 
a satisfactory understanding of properties of a particular model (especially on a 
critical line) one has to know distributions of various macroscopic (and, ideally, 
microscopic) quantities. For that reason in the present paper we study the limiting 
Gibbs states and the distributions of thermodynamic observables. 

One of the properties particular to disordered systems in statistical mechanics is 
the self-averaging of thermodynamic observables, introduced by Pastur and Figotin 
in the paper [IT]. There they also proved a general theorem concerning the self- 
averaging of thermodynamic observables for a wide class of models. By observables 
they meant quantities already averaged over the Gibbs distribution. For disordered 
systems involving a (realization of a) random field {hj,j G Z'^} the self-averaging is 
defined as follows. 

Definition (see [H]). A thermodynamic observable (Qn) is self-averaging, if 

lim (g^v) = Q (1) 

A'^oo 

exists and is the same for almost all realizations of the random field, where is the 
size of the system. 

The name self-averaging indicates that one does not have to average the ther- 
modynamic observable Qn over the distribution of the random field. Indeed the 
limiting distribution is concentrated at the average value, since Eq. ([1]) trivially 
implies 

lim {Qn) = EQ, 

N—*oo 

where E{-) denotes the average over the distribution of the random field. As a rule 
self- averaging observables are uniformly integrable, see [SI [15], hence, it is also true 
that 

lim E{Qn) = Q. 

Af— >oo 

From probabilistic point of view there are no fundamental differences between 
the thermodynamic randomness (described by the Gibbs distribution) and the ran- 
domness of the field {hj,j G Z"^}. Therefore it seems natural to get rid of the 
thermodynamic averages in the definition of self-averaging for observables like the 
magnetization. 
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Definition 1. A thermodynamic observable Qat is self-averaging, if 

lim Qn = Q, (2) 

N^oo 

exists and is the same for almost all realizations of the random field {hj,j G Z'^}, 
where the limit is understood in probability w.r.t. the thermodynamic randomness. 



There are thermodynamic observables which are not self-averaging on critical 
lines/points, having continuous (non-thermodynamic) distributions. For instance, 
it is widely known that the susceptibility 

1 ^ 

j,k=l 

is an observable of that kind. On the other hand, there are observables which 
distributions concentrate at a few (two or more) points. This fact motivated the 
authors of the paper [3J to introduce the notion of the conditional self-averaging. 

Definition 2 (see |3j). A thermodynamic observable Qn is conditionally self- 
averaging, if 

lim Qn — E{Qn\^n) = 0, in probability, (3) 

where E{-\C,n) are the conditional averages w.r.t. a sequence of functions of the 
random field {hj,j G Z'^} which obtain only a finite number of values, F, the same 
for all N. 

For an illustration of the notion of conditional self-averaging one can look at the 
random-field Curie- Weiss model, see [3]. In this model a conditionally self-averaging 
observable Qn is the magnetization J2f=i Sj, the sequence of functions ^n is the 
sign of the total random field 



N 




and E{Qn\^n) ~ ^n^*, where m* is the spontaneous magnetization. 

For a self-averaging observable Q n both thermodynamic (described by the Gibbs 
distribution) and non-thermodynamic (produced by the random field) fluctuations 
vanish as — ^ oo. It seems useful to introduce exponents which indicate how 
fast that happens. The exponent p, related to non-thermodynamic fluctuations, is 
defined by 

{QN-E{QN))=N-PrN, (4) 

as A^ —>■ oo, where the sequence of random variables vn converges to a random 
variable with a proper, non-degenerate distribution. The exponent r, indicating the 
magnitude of thermodynamic fluctuations, is defined by 

Qn - E{Qn) - N-^tn = N-^tN, (5) 
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as N ^ oo, where, again, the sequence of random variables t]\f converges to a random 
variable with a proper, non-degenerate distribution. The definitions of exponents p 
and r generalize straightforwardly to the case of conditional self-averaging. 

As a rule thermodynamic systems outside critical lines/points are collections of 
random variables {xj}jLi with short-range correlations. In this case one usually has 
self- averaging with the exponents p = t = ^. More precisely, 

1 ^ 

The exponents p and r are not fundamentally novel quantities. For most commonly 
used thermodynamic observables Qjy they are related in some way to the standard 
critical exponents. The values of exponents for the magnetization of the spherical 
model are calculated in this paper. 

Somewhat different terminology was used in the papers [H, [16]. There self- 
averaging with exponents p = r = | is called strong self-averaging, while self- 
averaging with exponents p G (|, 1) and r G (|, 1) is called weak self-averaging. 

Some general results on the behaviour of models under the influence of random 
fleld were obtained in the 70s and 80s by application of the renormalization-group 
ideas to the Ginzburg-Landau model, see [21 E]- In particular it was noticed that 
the random-fleld fluctuations dominate over the thermodynamic fluctuations as the 
critical point is approached. This observation suggests that the random-fleld fluctu- 
ations also dominate on the low-temperature critical line, and hence one should have 
T > p there. This is exactly what happens with the fluctuations of the magnetiza- 
tion of the spherical model, and we will see in Section 6 that in this case P = | ~~ f 
and r = I - i. 

The rest of the paper is organized as follows. Section 2 contains the exact 
deflnition of the spherical model, the random fleld and the boundary conditions. 
It also contains some well known technical results for the use in the later sections. 
Section 3 summarizes the main results of the paper. In Section 4 we calculate the 
free energy of the spherical model as an illustration of the application of saddle- 
point method in the low-temperature region. In Section 5 we describe in details the 
properties of the spherical model (the random fleld {xj,j G Z'^}) in the inflnite- 
volume limit. In Section 6 we provide an analogous detailed description for the 
magnetization of the spherical model. The results of Sections 5 and 6 in the absence 
of the boundary fleld are re-derived in Section 7. The results of the paper are 
discussed in Section 8. 

2 The model and useful facts. 

The spherical model describes a collection of random variables {xj,j G Z'^} placed 
at sites of an integer c/- dimensional lattice, Z'^. Every site j G Z'^ is specifled by its 
d integer coordinates (ji, j2, ■■■,jd)- 
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To define the distribution of random variables at all sites of the lattice, we first 
specify the joint distribution for the random variables in a finite rectangle 



Vn^ljeZ^-.lK <n,iy^ 1,2, d} 



containing N = n'^ sites, and then pass to the limit n —>■ oo. To avoid unnecessary- 
complications we impose periodic boundary conditions in dimensions 2,3, ... ,d. 
Thus the boundary of the rectangle Vn is the set 



Bn = {j e Vn : ji = l,n}. 



The Hamiltonian. 



The random variables located in the rectangle Vn interact with the boundary field, 
the external random field, and each other via the Hamiltonian 

Hn — J ^ ^ Tji^XjXk ^ ^ ^j'^j ^ ^ ^ ■^j-i 

where J > 0, Tjk arc the elements of the nearest-neighbour interaction matrix, 
{hj, j G Z'^} is a fixed realization of the external random field, and b is the boundary 
field. 

The interaction matrix. 



The elements of the interaction matrix T are given by 

1^=1 ie{i,2,...,d}\i' 



where 



S{ji,ki 



1, \{ji = ki, 
0, liji^h. 



is the Kronecker delta. 

The coefficients J'^^\j\, ki) are the elements of the n x n tri-diagonal matrix 



/o 1 

1 

2 











5 oy 
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The coefficients J'^''\jt,, /ov), for u = 2,3, . . . , d, are the elements of the matrices J^'^^ 
which have extra | at the upper right and lower left corners (due to the periodic 
boundary conditions) 



/o I 

i i 











1 

2 



The eigenvalues of the matrix J^^^ are given by 

A; = cos , l — l,2,...,n. 

n + 1 

The corresponding orthonormal (that is, orthogonal and normalised) eigenvectors 
are given by 



(0 



2 . TrZm 
sm 



n + 1 n + 1 



, Z = l,2,...,n. 



m=l 



The eigenvalues and orthonormal eigenvectors of the matrices J'^'^\ for u = 2,2>, . . . ,d, 
are given by 



27r(Z - 1) 

A; = COS — ^ I = 1,2, . . . ,n, 



n 



and 



u 



(0 



(0 



COS 



27r(Z- l)(m- 1) n 
n 4 



, / = l,2,...,n. 



m=l 



Finally, the eigenvalues of the interaction matrix T are the sums of the eigenvalues 
of the matrices J^''^ 

d 

Afe = Afc^ + ^ Afc^, /c = {h, k2,..., kd) e K- 

The corresponding orthonormal eigenvectors are the products of the eigenvectors of 
the matrices J^'^^ 



w 



I i/=2 



, k = {ki,k2,...,kd)eVn. 



(6) 



' jeVn 



The external random field. 
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We assume that the coefficients {hj,j G Z'^} are a fixed reahzation of indepen- 
dent normal random variables {hj,j G Z'^} with zero mean and variance h"^. The 
assumptions of independence and normal distribution are made to avoid unnecessary 
complications. The behavior of the model is very different if the random variables 
{hj,j G Z'^} have, say, Cauchy distribution, or, if the random variables have strong 
negative correlations severely suppressing fluctuations of sums like J2j£V„ ^j- Nev- 
ertheless, we restrict our attention to the technically convenient case of independent 
normal random variables where the fluctuations are neither abnormally large, nor 
abnormally small. 

The Gibbs distribution. 

The distribution of the thermodynamic random variables {xj,j G Vn} is specifled 
by the usual Gibbs density 

p{{Xj,J G Vn}) = -Fj— , 

with respect to the spherical "a priori" measure 

Hn{dx) = E ~ ^ n ^^j- 
The normalization factor (partition function) 0„ is given by 

/OO nOO 
... e-f'''-fin{dx). (7) 
-OO J — OO 

Useful estimates. 



Equations (l8l)- f[T2|) below state well known results which are used throughout the 
paper. A routine analysis of the singularity at lui = a;2 = . . . = cu^ = shows that 
the function 



1 



/ w \ J. i 



1^ duj^ 



< OO 



at z = d if d > 2m. 

Let 7 G [0, 2), ^ > 0, and Zn = Amax + C"-"''') then we have as n — »• oo 



AT ^ 

k&V-a 



^ W^r\^n) - ^/\wt\zn) + o [exp (-ni-^/^c(C) 



2n 



where 



/\Wt\zn) ^ Wtl{Zn - 1) + Wtl{Zn + 1) - 2W!r\znl 

and c(C) is strictly positive and increasing for C > 0. 



(9) 
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If 7 = 2, C > 0, and > 4, then 
^ E 't X = ^i'^(^) - ^^^i'^(^) - CW^i'^(^) + o(n-2), (10) 

as n oo, where the prime indicates that the summation does not involve k = 
(1,1,. ..,1). 

li d > 2m, and C ^ 0; then 

i E ^,1-2 = Wt\d) + 0(1), (11) 

as n ^ oo. Approximation of sums of the type ([9]), ( |TOl) by integrals was analysed 
in [H [7] . For an outline of a method particularly suited for the above sums see [12] . 
If m > and d > 2m, then 

. exp(^EtiX.a;.) du. r(d/2 - m) /f .r"'/^ 



as X]l5=i^^ C)0- Foi' derivation of the above asymptotic formula in the case 
m = 1 see, e.g., [9]. The method used in [9] can be also apphed in the case m > 0. 

Finally, a direct numerical computation of the multiple integrals wj[^\z) is an 
awkward task. Fortunately, for m > and z > d, it is reduced to the following 
integral of the Bessel function /o(a;): 

/ \ 1 /"oo 



3 The main results. 

Usually thermodynamic properties are derived in the limit of an infinitely large 
lattice. In our case the results are most conveniently formulated in a continuum limit. 
We choose to use the version of continuum limit where the limiting configurations 
are random functions defined on the (i-dimensional rectangle [0, 1]'^: 

Ml)}^em^ = {a;(7i,72,---,7d)}7i,72,...,7de[o,i]- 

For any 7 G [0, 1]*^ the random variable ^(7) is defined as the following limit in 
distribution 

where [y], is the integer part of y. 

Thermodynamic random variables ^(7) and x{6) are limits of the random se- 
quences X(^[^-^n],['Y2n],...,l-ydn]) ^ud X(^[s-^n],[52n],...,[5an]) Separated by a distance of order n. 
Hence, in the continuum limit the random variables x{'y) and x{6) with 7 7^ 5 are 
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independent due to the exponential/power-law decay of thermodynamic correlations 
in the high/low temperature region. 

Unless explicitly stated otherwise, in this paper we consider dimensions > 5 and 
inverse temperatures j3 > /3c, where (3~^ is the critical temperature of the spherical 
model in external random field, see the paper by Pastur |13j . 

Denote v^(i,i,...,i) the projection of the external random field on the eigenvector 



w 



(1,1,.. .,1) 



corresponding to the maximal eigenvalue of the interaction matrix 



<^(i,i,...,i) 



n 



(n + 1) , 



} ^ sm r III. 



n + 1 



Recall that the external field {hi, I G Z"^} is a realizations of the random field 
{hi, I E Z'^} of independent normal random variables with Ehi = 0, Ehf = h?, 
for any / G Z*^. Everywhere below we will use the notation ^{a, b"^) to denote 
thermodynamic normal random variables with mean a and variance 6^, which are 
independent from the external random field. The symbol q will be used to denote 
a realization of a non-thermodynamic normal random variable q. The value of q is 
fixed once we fix a realization of the random field, and q is always independent of 
M{a,b^). 

The main results of the paper can be stated as follows. 

1. In the absence of the boundary field, 6 = 0, the random variables x{'y) have 
normal distributions with the expected values 



(x(7)) = sgn 
and the variances 



V5(i,i,...,i) 



sin(7r7i 



+ ^7' 



where are independent realizations of zero-mean normal random variables 
with the common variance 




2. For a fixed realization of the external random field, the law of large numbers 
is valid for the normalized sums 

as n —>■ oo. The convergence to the limiting value can be summarized by the 
following asymptotic formula 



m„ ~ sgn 





2 


( f3AW^'\d) 


_<^(l,l,...,l) 







+ n^-'^'^qn+ 
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n 



-d/4 



Vlv^{i,i,...,i) 



where g„ is a reahzation of a zero-mean normal random variable with the 
variance 

„ Ttt^ - 69 / /i \ ^ 



37r6 



2J 



Hence, the magnetization m„ is (only) conditionally self- averaging with the 
exponents P = \ — 2 ^^"^ ^ ~ \- 

3. For 6 7^ the random variables x(7) have normal distributions with expected 
values 

6cosh[(l-27i)v^ 

(x(7)) 



J cosh v(o 



and variances 



where Co is a solution of Eq. (fT6l) . and are independent realizations of zero- 
mean normal random variables with the common variance 



h 



4. For b ^ 0, the law of large numbers is valid for the normalized sums 

'^n = ly X! -^i' 

as n — > 00. The convergence to the limiting value can be summarized by the 
following asymptotic formula 

m„ ~ + n^-"\. + n'-'/W„ fo. fl - 

where g„ is a realization of a zero-mean normal random variable with the vari- 
ance (1301) . Hence, the magnetization m„ is self-averaging with the exponents 



tanh ^/(^ 



P = ^ - i and r 



1 _ 1 

2 d- 



4 The free energy. 

The calculation of free energy, expected values and correlation functions for the 
spherical models is reduced, in a routine fashion, to calculation of the large-n asymp- 
totics of an integral. In this section we find the large-n asymptotics for the free 
energy 

fn = — TT-J In ©n- 
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A particular attention will be paid to 0(n^^) asymptotics of /„, which, as it turns 
out, determines thermodynamic properties of the model below the critical temper- 
ature. 

The introduction of new integration variables {yj}j£Vn Eq. ([7]) via the orthog- 
onal transformation 

= Vfc, j e K, 

k£V„ 

where the eigenvectors {'Wj''^}j^v„ are given by Eq. ([H]), diagonalises the interaction 
matrix. Therefore, we obtain the following formula for the partition function 



where 



k€Vn k£V„ k€Vn 



^k 



Since the vectors {Wj''^}j^v„, k & Vn are orthonormal, the random variables ipj^ = 
J2jev„ hjwf\ are independent normal random variables with zero mean and vari- 
ance /i^. Therefore, we can treat the coefficients (pk-, k G Vn a.s realizations of 
independent normal random variables. 

A direct calculation of the coefficients a^, k = {ki, k2, ■ ■ ■ , kd) & Vn (using only 
the formula for the sum of a geometric series) yields 



(fc) 



and Qfc = 



w 



(k) 



n+l 



5{k2,l)...5{kdA) X 



sm 



irki 
n + l' 
0, 



if ki is odd, 
if ki is even. 



The integral representation for the delta function 




N 



/ ds exp 

21X1 J-ioo 



in the "a priori" measure allows one to perform integration over the variables yj, 
j G Vn- However, we can switch the order of integration over the variables yj, j & Vn 
and s only after a shift of the integration contour for s to the right. The shift should 
assure that the real part of the quadratic form involving the variables yj, j & Vn is 
negatively defined. The switching of integration order, integration over yj, j G Vn, 
and the introduction of a new integration variable z via s = (3Jz yields 



13J 
27ri 



TT 



N/2 



+ioo+c 



joo+c 



dz exp [N/S^niz]] 



(13) 
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where 



^^^^ fcey„ 4-^^^ fcGVn ^ - ^fc 



\2 



and c > (i is the shift of the integration contour mentioned above. 

The large-n asymptotics of the integral f|T3|) can be found using the saddle-point 
method. The saddle point of the integrand is a solution of the equation 

For any z > d, as n ^ oo, the sequence of the derivatives $^(-2) converges, with 
probability 1, to 

^'{z) = j-^w^'\z)-^w:t\z), 

where the functions wl[^\z) are defined in Eq. ([8]). The function $'(2) increases 
monotonically with z on [d, 00), and the location of its zeroes depends on the dimen- 
sion d of the lattice. Namely, if d < 4, then the function $'(2;) has exactly one zero 
on the interval [d, 00) at a point z* > d, for any (3 > 0. If > 5 and the variance of 
the external field, h"^, is sufficiently small, then there exists a critical value 

1 



2^1-(A)Vf(c^) 



1 - f , 

of the parameter P, see pjj. If /5 G (0,/5c) (the high-temperature regime), then the 
function $'(2;) still has exactly one zero on the interval [d, 00) at a point z* > d. 
While if (3 > (3c (the low-temperature regime), then the function ^'{z) is strictly 
positive on the interval [d, 00). 

The application of the saddle-point method for the integral (fT3l) is fairly straight- 
forward when the saddle point z* is greater than d, see [5]. Therefore, in this paper 
we consider only the low-temperature regime and d > 5. When (3 > (3c-, the function 
^n{,z) still attains its minimum on the interval (Amax, 00) at a point 2;* > Amax, 
where Amax = d — 1 + cos is the maximum eigenvalue of the interaction matrix 
T. However, the sequence of saddle points z* approaches the branch point of the 
integrand at z = Amax, and the application of the saddle-point method becomes a 
bit more tricky. 

To be able to apply the saddle-point method we have to find a change of variables 
z = Amax + C'^"''') such that the sequence of rescaled saddle-points = (2* — Amax)^''' 
converges to a positive limit > as n ^ 00. Then, the application of the saddle- 
point method for the integral over C, becomes straightforward again. Note that the 
above search for a proper change of variables has an important physical meaning — 
rC^I'^ l\fC^ is the correlation length of the model. 

In order to find the proper value of 7 we have to analyse the sums in Eq. ([1] 
The large-n asymptotics of the sum 
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when z = Amax + C^ and C > 0, follows from Eqs. ([9]) and (ITO!) . Namely, as n — >• oo, 



To find the large-n asymptotics of the sum 



when z = X^ax + C^"'''? we have to use the law of large numbers. First, we take out 
the term corresponding to/c = (l,l,...,l) and rearrange the sum as follows 



For C > 0, Eqs. ([9]) and <^ yield as n oo 

Let {ij,n\^=i^=i be a triangular array of independent random variables with zero 
expected values. The condition 

n 

E E\S,j^nY 0) for some s G (1, 2], 
i=i 

as n ^ oo, is sufficient for the validity of the law of large numbers 

n 

XI ^hn — ' 0, in probabihty, 
see, e.g., [B]. Therefore Eqs. dH]), dH]), and (fTTl) imply 

T7 'tt -V^:; 0' in probability, 

as oo, if (i > 4, and C ^ 0. Summarizing the above we obtain 

2 

S2(A,„ax + C^-^) = + h:'Wf\d) + 0(1), 



as n — > oo. 
The sum 
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with z = Amax + is a realisation of a normal random variable with zero mean 

and the variance 

It is possible to find a relatively simple expression for the variance 



2h' ^ (l + (-l)'^+l)^sin^„^^ 



\ n+l ^ 71+1 



First, note the identity, see |12 



where 



x{z) = l + z-d + {z - d){2 + z - d). 
On differentiating Eq. f[T^ over z three times we obtain 

m (z- Afc)4 ~ n'^+i (z-ci)(2 + z-d) [[{z-d){2 + z- d)]^/^ x^+\z) + 1 
4(n + 2(r2 + l)a;("+i)(^) /n(n + 2) 1 



+ 



Hence, if 7 G (0, 2), then 



^d+l-57/2 (2()5/2' 



as 77, — >• 00, while if 7 = 2, then Amax + C"- ^ ~ — ^Tc'^n ^ + aiid ct^(C) 

n^~'^h^t{C, — ^TT^), where 

1 / 1 / — 1 2C — 3 

^(0 = 7 TTTw^tanh JiC wf= + 



C I (20^/^ V 4 cosh^ M cosh^ a/IC 



The function t{C,) (and similar functions below) has only a removable singularity at 
C = 0, and the analytic continuation is to be used for negative values of C,. Thus 



S3 (Amax + C'^ ^) ~ O (n 



57/4-(c(+l)/2 



does not produce a non- vanishing contribution to the saddle-point equation if 7 < 2. 
It is also possible to obtain a simple formula for the sum 
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by differentiating Eq. (fT5|) over z. The differentiation yields 



n 



n+l| 



+ 1 x^f^;) 



+ 



On replacing z by Amax + C'^ ^ we obtain 

4^-1+7/2 



-C 



2tanhyi(C" 



+ 



cosh%/i(C-V^ 



if 7 e (0,2) 



if 7 



Thus the sum T,/^{z) is dominant among the four sums S;(2;), / = 1, 2, 3, 4 (if 6 7^ 0), in 
the sense that it is 24(2) that controls the location of the saddle point in the low- 
temperature region. Indeed, the sum T,i{z) produces a non- vanishing contribution to 
the saddle-point equation already in the scale z = Amax + C/'^^- Moreover, the extra 
contribution produced by Ti^i^z) prevents the rescaled saddle-point (* approaching 
the branch-point at C = 0, where the remaining sums could, potentially, yield non- 
vanishing contributions to the saddle-point equation. 

On introduction of the new integration variable ( in Eq. ( |T3|) via z = Amax + C^~^ 
we obtain 



,-1 



/3J 



N/2 



ioo+Co 



dC exp Np<i>n{\me.^ + Cn- 



The saddle-point of the integrand is (* = 2(q + ivr^, where Co is a solution of the 
equation 



h 



2(3J 



b \ ^ f tanh y/Q^ 



2J 



+ 



'^J ) \ VCo cosh^ 



Application of the saddle-point method yields 



(16) 



' fn 



pn"- 



— In — + $(rf) + n- Vi + n-'(j)2 (2Co) + o{n 



as n — i> 00, where 



... j In z - ^ COS cj;, 
^ AL,(t;)-^AW^i''(rf) + ^' 



M 27r' 



02(0 = {j 



8J 



J' 



62 

J 



2Ctanhv/|C 
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The function ^{z) determines the thermodynamics of the model in the high- 
temperature region. The term 0i appears because of to the lack of periodicity in 
one of the dimensions. The function 02(C) is responsible for the thermodynamic 
properties of the model on the low-temperature critical line. 



5 Individual distributions. 

To find the individual distributions of the random variables {xj,j G V^} we calculate 
the corresponding characteristic functions 

K,j{t) — {exp{itXj)). 

The saddle-point method described in the previous section yields the following large- 
n asymptotics 



Kj{t) ~ exp 



w 



E 



(17) 



Therefore, for large values of n, the individual distributions of the random variables 
{xj,j G Vn} are nearly normal with mean values 



keVn - ^k 



(18) 



and variances 



2pJ 



fceVn n ^k 



On substitution z*^ — (*n ^ one obtains 



2/3J 



/TT /"TT 
-TT J — TT 



dui 



d - T,t=i cosu,, f-Ji 27r 



n 



as n ^ oo. Thus, in the low-temperature region and in the presence of the boundary 
conditions, b ^ 0, the variances of the thermodynamic random variables Xj are not 
affected by the random field {hi, I E Z^}. As ji increases. 



COs(jiCJi' 



cos(jii 

J-TT J-TZ d — J2t=l 



-ji duJi T{d/2 - 1) 



COS uj^ 27r 27r°'/2j^ 



d-2 



0. 



Hence, only random variables near the boundary have variances noticeably different 
from the bulk value 



'bulk 



Wi'\d). 



2PJ 
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The first half of the sum in Eq. ([18 



— r, T 2^ 



describes the shift in the expected value of Xj due to the external random field. It 
is a realization of a normal random variable with zero mean and variance 



2 / „.,ik) 



h \ ^ w 



3 



As n — i> oo the variance V^' tends to 



1 — cos(jiCJi) duji 



For d > A we have 



l-n'" {d - ELi cosu;,)2 27f 47r'^/2jf-^ 
as ji — > oo. Hence, the variance V^^ also approaches its bulk value 



as we move away from the boundary. 

The second half of the sum in Eq. (fTSjl . 



b „ akw) 



"3 



be ^ ^ ^ 



k&Vn 



is the shift in the expected value of the thermodynamic random variables Xj due to 
the influence of the boundary conditions. An application of the "contour summa- 
tion" technique, see [12], yields the following simple formula 

The large-n limit of fi^ depends on the location of the node j = {ji,j2, ■ ■ ■ ,jd)- 
Assuming ji ~ •jin as n — > oo, we obtain (recall that = Amax + C*^"^ ~ d+2(on~^ 
in the low-temperature region, see Eq. ( |T6l) ) 



6 cosh [(l-27i)v/Co] 
hm /i'"^ = ^ — - — 7= ^ = fj. hi). 



J cosh 
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The characteristic function of an arbitrary pair {xj,xi) is given by 

ts 



K,j^i{t, s) = {exp{itXj + isxi)) ~ Kj{t)Ki{s) exp 



(k) {ky 



as n ^ oo. Hence, for large values of n, the joint distribution of Xj and xi is nearly 
normal with the covariance 



COv{Xj,Xi) 



1 



2PJ 



(k) (k) 
W j vo l 



fcey„ ^k 



E 



Since z* 



Amax + C*'^ ^) "we have (ignoring thin layers near the boundaries) 



lim cov {xj,xi) 



1 



2/3 J J-n 



exp 



COStUi, 
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Thus, the covariance cov(xj, xi) shows the usual, for the critical line of the ordinary 
spherical model, power-law decay with the distance r'ji = J2f,=iUi' ~ ^uY between 
the nodes j and I. Indeed, using Eq. (fT2l) we obtain 



cav{xj,xi) 



T{d/2 - 1) 



(20) 



if 1 <^ Tj^i <ti n. 

Summarizing, we conclude that the structure of random variables {xj^j G ¥„} 
is fairly simple. Ignoring thin layers near boundaries, we have in the limit n — > oo 



(7i),<ik), 



where qj is a realization of a (non-thermodynamic) normal random variable with 
zero mean and the variance V^^^ii ^"^^ Mj{a, 6^) is a thermodynamic normal random 
variable with the mean a and the variance 6^, see Fig. 1. 

In the presence of the boundary conditions, apart from the global influence 
through the saddle point C,* , the external random field {hj^j G Vn} produces 
only additive contributions (random shifts) qj to the thermodynamic random vari- 
ables G Vn}- The properties of the (non-thermodynamic) random variables 
{^jij ^ generating the shifts are fairly interesting. At the critical temperature 
(3c the random field {Qjij G Vn} undergoes a transition into a phase with long-range 
correlations, see Fig. 2. 

Indeed the covariances of the random variables {^j, j G Vn} are given by 



cov(qr.,qr;) 




Passing to the limit n 



Jimcov(q,,9,)= — 



oo we obtain 



E 

keVn 



exp 



(fc) (fc) 
w(j 'vol 



COSU, 



n 

i/=i 



~2^ 
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Figure 1: A line of a 5-D realization of thermodynamic random variables {Xj,j G 
Vn} (discs) for (3 > jSc- The picture also contains the corresponding realization of 
the random field {qjij G Vn} (circles) driving the random variables at very low 
temperatures. 




Figure 2: A line of a 5-D realization of the random field {q^, J G Ki} for j3 > j3c- In 
comparison with independent random variables the field {q^, J G Ki} has a substan- 
tial inertia — positive/negative values tend to be surrounded by positive/negative 
values. 
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If /? < f3c, then z* > d and the above integral decays exponentially with the distance 
Tj^i between the nodes j and /. If P > Pc, then z* = d and Eq. (fT2l) yields the 
power-law decay 



h^T{d/2 - 2) 



Note that the correlations of the random field {qj,j € Vn} decay noticeably slower 
than the correlations of the thermodynamic random variables, see Eq. (1201) . This 
slow decay of the covariances is the reason for the dominance of the random-field 
fluctuations over the thermodynamic fluctuations. 



6 Macroscopic observables. 



Our aim in this section is to establish the law of large numbers for the normalized 
sums (magnetization) 

= ly II ^i' (21) 

and to study fluctuations (the central limit theorem) of these sums around the 
limiting value. The corresponding characteristic functions are given by 

The large-n asymptotics of K„(t) is calculated using the technique of the previous 
section. The saddle-point method yields 

where (see Eq. ([6])) 



for k = {ki, k2,..., kd) G K- 

Thus, for large values of n, the distribution of the magnetization (|2T1) is approx- 
imately normal with the mean value 

The sum 

1 ^ ^kVk n~"y.(,,,i„„,i)[l -(-!)'] sin^ 



^JN,i^„z*-Xk ~ 4J V^i^r-f-z^-ci-cos^l-cos^' 
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is the shift in the expected value of the magnetization (pT]) caused by the external 
random field. It is a realization of a normal random variable with zero mean and 
the variance 



4J2(n + 



sm 



2 nl 
n+1 



COS ■ 



{l + z*-d 

On calculating the sum over I, see [12], we obtain 



cos 



nl \ 
n+1 J 



2 • 



s„ 



n 2x"+^«)+a:"(z:)-x(4)-2 ^ (n + l)x"+^(4) 
2 + 1) {x{zl) - x-\zl)) {x^+\zl) + 1)' 



(24) 

In the low-temperature region we have = A^ax + C*^ ^ ~ + 2C,Qn ^, see Eq. 
(flGl) . therefore 



h \ n 



2 



2J; 4Co2 



^ _ 3 tanh VCq _^ 



(25) 



as n — >• cxo. 
The sum 



■ 2cosh^v^y 

Qfc^fc ^ b 2x{zl) x"«) - 1 
2'^^ < - Afc ~ Jnx(z*) - la;"+i(<) + 1' 

is the shift in the expected value of the magnetization (!2T!) caused by the boundary 
conditions. On substitution z^ = Amax + (*n~'^ one obtains 



E 



b ^ afc?7fc b tanh v^Co 



(26) 



as n ^ cx). 

Let's now look at the variance of the magnetization. According to Eq. (!22|) it is 
given by 



1 



E 



n 



-d-l 



2PJN^- zl - Xt 2pj(n + 1) 1 + < - d 



1 = 1 ^ 



1 



-11 



sm 



2 ttI 

n.+l 



COS ■ 



n+1 ( 1 — COS 



n+1 



2 • 



The remaining sum over / can be calculated exactly, and we obtain the following 
expression for the variance 



n 



'd-i 



a 



2pJ{z*^ - d) 



n 



2x{zl) {x-{zl) - 1) 



ixizl) - 1) ix^'+^iz*) + 1) 



(27) 



On substitution of 2;* = Amax + C*'^ ^ for the saddle-point one obtains 

o ( tanh ^/Co \ 



(7 ~ 



4/?JCo 



(28) 
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Figure 3: The infinite-lattice magnetization m = linin^oo'^n as a function of tfie 
normalized boundary field b/J, for d = 5, (3 J = 2, and h/J= 0.5. The left/right 

2^2 

limits at 6 = are given by =F 

vr 




as n — > oo. 

Summarizing the above we obtain the following expression for the magnetization 
m„ ~ 4^^^ + n-'^/^g^ + n^-^^^Af. (o,-^(l- '-^^] ] , (29) 



where g„ is a realization of a zero-mean normal random variable with the variance 
fhY 1 3tanhv/<o^ 1 \ (3^^ 



V2j; ACl V 2v/Co 2cosh2 

and Afnifi, v"^) is a thermodynamic normal random variable with mean /x and variance 
f ^. Therefore, the magnetization of the spherical model is self-averaging (for & 7^ 0) 
with the exponents P = | ~ | and r = | — ^- The limiting magnetization m = 
limrj^no function of the boundary field h is shown on Fig. 3. 



7 The distributions for zero boundary field. 

As it is clear from previous sections a non-zero boundary field dominates over the 
(zero-mean) random field in the low-temperature regime. Therefore in this section 
we consider the case of zero boundary field. 

If 6 = 0, then the saddle-point equation for the integral (fT3|) is given by 

*«'^'-^-^S7^-ji^E(^)^^o. (31) 



22 



Again, the saddle-point 2* drifts towards the branch-point of the integrand in the 
scale where the terms corresponding to k = (1,1,...,1) produce a non- vanishing 
contribution to $^(-2). As is obvious from Eq. (^11) . that happens in the scale 
z = ^(1,1,. ..,1) + C'^ The distance from the saddle point 2;* = A(i_i^...^i) + (*n~'^^'^ 
to the eigenvalues Xk with A; 7^ (1, 1, . . . , 1) is at least of the order 0(n~^). Therefore 
there are no additional non-vanishing contribution to the saddle-point equation from 
those eigenvalues. 

In the scale z = A(i + (n^'^^'^ we obtain the following saddle-point equation 
in the limit n ^ 00 



1 



1 



h 



(2), 



2PJ 



1 

4J2 



0. 



The positive solution of the above equation is given by 



2JW 1 



The location of the saddle-point z*, as n ^ 00, is given by 



zl^d-\ 



TT 



COS ■ 



n 



Evaluation of the characteristic function (1171) at 2;* shows that the thermody- 
namic variables Xj have normal distributions with the expected values 



1 V5(i,i,...,i)«^, 



(1,1,.. .,1) 



2J 0* - A(i,i„„,i) 



- E 

2J ^ 



(fc) 



fcev„\{i,i,...,i) ^« 



and variances 



2/?J 



E 



w 



(fc)y 



Assuming that j = (ji, J2, ■ ■ ■ ,jd), and that for k = 1, 2, 
with 7fc G (0, 1), we obtain 



d we have jk ~ 'Jk'n 



lim Hj = sgn v?(i,i,...,i) sin(7r7 



(32) 



where 7 = (71, 72, ... , 7^), and are realizations of independent zero-mean normal 
random variables with the variance V^^Yk given by Eq. (1191) . An important feature 
of Eq. fl32|) is the term sgn(y9(i^i^...^i) common to all expected values fij. This term 
is the reason for the absence of conventional self-averaging for the normalized sums 
(magnetization) 
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On substitution of the saddle point 2;* = A(i^i^.,,^i) + (*n~'^/'^ in Eq. (!22l) we see 
that, as n — >■ 00, the distribution of the magnetization (133|) is asymptotically normal 
with the expected value 



E 



+ 



E 



~ sgn [V5(i,i,...,i) 
fkVk 



\ 



. /3 /3cJ 



^"^^ fcey„\(l,l,...,l) ^ri ^k 



On subtracting the contribution of the maximum eigenvalue from Eq. fl2^ one finds 
that the remaining sum over A; is a realization of a normal random variable with 
zero mean and the variance 



i-d 



as — 00. 

On substitution of the saddle point = A(i^i^...^i) + (*n~'^/'^ in Eq. fl271) we find 
that the thermodynamic variance of the normalized sums fl33|) is given by 



^2 



1 



vr2|y.(i,i,...,i)|/3n'^/2\ 



as n — >• 00. 

Summarizing the above we obtain the following expression for the magnetization 



nir. 



+ 



sgn 



n 



V5(i,i,...,i) 

-d/4 



1 - 



/3 / /3,J 



+ n 



2~d/2 



qn+ 



V 



TT 



1 - 



(34) 



where g„ is a realization of a zero-mean normal random variable with the variance 

„77r2-69 ( hV 



37r6 



2J, 



and A/'„(0, f ^) is a zero- mean thermodynamic normal random variable with variance 
f ^. Thus, in the absence of the boundary field, the magnetization of the spherical 
model is conditionally self- averaging with the exponents P = \ — 2 ^"^^ ^ ~ \- 



8 Discussion and concluding remarks. 

It was shown in the paper [3] that there are problems with almost sure convergence 
of Gibbs states for the random-field Curie- Weiss model in the infinite- volume limit. 
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In fact, below the critical temperature, the limits of thermodynamic averages {sj)N 
do not exist, almost surely, as the volume N tends to infinity. A possible solu- 
tion of the convergence problem was also proposed: it is necessary to consider the 
limits of distributions of {sj)]\f, which, after some minor technical efforts, lead to 
correctly defined random infinite-volume Gibbs states. The same problem exists in 
the spherical model, and, most likely, in such often considered models as the Ising 
model and 0{n) models. Namely, for 6 = 0, lim„^oo(3;j)n does not exist almost 
surely, although it exists in distribution. The results of the present paper show that 
switching on a homogeneous boundary field rectifies the problem with almost sure 
convergence. Namely, for b ^ 0, \imn^oo{xj)n exist almost surely, which (together 
with convergence of higher correlation functions) means that the corresponding limit 
Gibbs state exists for almost all realizations of the random field {hj, j G Z'^}. 

At approximately the same time, Newman and Stein [TDj pointed out that the 
absence of convergence of local thermodynamic averages, like {sj)^, is a natural 
occurrence in many disordered systems. They call this phenomenon the chaotic size 
dependence. Somewhat later, Newman and Stein also proposed their own solution 
of the problem with infinite- volume Gibbs states. Instead of looking at distributions 
of local averages like (sj)jv they choose to look at the empirical distributions 

J^^iy) ^ N-'i^{k G {1, 2, . . . , AT} : {s,)^ < y}, 

for a fixed realization of randomness. Assuming ergodicity we have 

lim J^N{y) = lim Pr [{sj)^ < y] , 

hence, both constructions provide the same result: a random infinite-volume Gibbs 
state. Newman and Stein call the random Gibbs state the metastate. 

The authors of the paper p] investigated self-averaging using the ideas of renor- 
malization group theory. They concluded that there are universality classes of mod- 
els within which a particular non-self-averaging thermodynamic observable has the 
same distribution in the thermodynamic limit. The results of the present paper 
indicate that the conclusion of the paper [T] looks plausible, at least for the mag- 
netization. Indeed, according to Eq. fl5^ the magnetization of the spherical model 
obtains the values ±m* with probability |, where m* is the spontaneous magneti- 
zation. The magnetization of the Curie- Weiss model and, most likely, of disordered 
finite-dimensional Ising models has the same distribution, see [3]. One can also guess 
that the magnetization of various disordered 0{n) models is uniformly distributed 
over an n-dimensional sphere. On the other hand, we also saw that the distribution 
of the magnetization is highly- sensitive to symmetry-breaking pertnThaiions. Indeed 
an arbitrarily weak symmetry-breaking boundary field restores self-averaging, that 
is, changes a non-degenerate distribution to a degenerate one. Although that fact 
rather goes along with than contrary to the lines of renormalization group argument. 

The susceptibility of the spherical model 

= /5n'^t-Var(m„), 
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can be easily found from Eqs. (!29|) and flM|) . If 6 7^ 0, then (when properly normal- 
ized) the susceptibility is self-averaging 



The susceptibility of various 3D disordered models was studied intensively using 
Monte-Carlo simulations since mid-90s, see, e.g. p3l [T6]. The histograms obtained 
in [m |16] suggest that the distribution of the susceptibility is not normal, positively 
skewed, and has heavy tails. The distribution of the susceptibility given by Eq. fl35l) 
has the same properties, and thus, to some extend, explains the results of Monte- 
Carlo simulations. It has been suggested in the paper that the distribution of 
susceptibility should be the same within universality classes. Since Eq. ( l35l) is the 
asymptotics of Eq. (1271) at the pole z = X^ax it is not unreasonable to expect the 
universality of the distribution of Xn for a certain class of models. Although it 
is tempting to speculate that 0{n) models might belong to the universality class, 
nevertheless, the results of the present paper do not indicate neither how wide the 
universality class is, nor which models possibly belong to this class. 

In conclusion, various disordered models have been intensively studied recently 
either numerically or using various heuristic approaches like, for instance, the renor- 
malization group. The present paper derives explicitly distributions of various ther- 
modynamic quantities within a non-trivial disordered finite-dimensional model — 
the spherical model in a random field. The author hopes that the paper is helpful for 
understanding the conclusions of heuristic theories, and for interpreting the results 
of Monte-Carlo simulations. 
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